WASSERSTEIN DECAY OF ONE DIMENSIONAL JUMP-DIFFUSIONS 

BERTRAND CLOEZ 



Abstract. We are interested by a one dimensional Markov process wiiiclr moves following 
a diffusion for some random time and then jumps. It can represent some natural phenomena 
like size of cell or data transmission over the Internet. The paper begin with some results about 
Lipschitz contraction of semigroup. Our approach is connected with the notion of curvature 
introduced by Ollivier and Joulin. Our main results for jump-diffusions are quantitative estimates 
in Wasserstein distance, when the jump times depend of the space motion. We use different 
techniques which are a particular Feynmann-Kac interpretation and a non coalescent coupling. 
Several examples and applications are developed, including explicit formulas for the equilibrium, 
application to branching measure-valued processes and integrals of compound Poisson process 
with respect to a Brownian motion. 
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1. Introduction AND MAIN RESULT 

1.1. Introduction and related work. We are interested by a process that models a fragment 
which remains stable and moves for some random time and then jumps. It represents some nat- 
ural phenomena that can be observed at a great variety of scales. To give just few examples, let 
us simply mention the studies of stellar fragments in astrophysics, fractures and earthquakes in 
geophysics, degradation of large polymer chains in chemistry, size of cell in biology, fragmen- 
tation of a hard drive in computer science. More precisely. This process X = {Xt)t>o has E, a 
subinterval of M, as state space and its infinitesimal generator is given, for any smooth enough 
function / : £' i-> M, by 

(1) Vx > 0, Cfix) = aHx)f"{x) + g{x)f{x) + r{x) f {F{x, 6)) - f{x)de^ 

where a, r, g and x i— )■ F{x, ■) are C°°, 9 i— >■ F{-, 9) is measurable, r and a are non negative. 
Betweens the jumps, this process evolves like a diffusion {Yt)t>Q which satisfies the following 
stochastic differential equation : 

(2) dYs = g{Ys)ds + V2a{Ys)dBs 

where {Bt)t>o is a standard Brownian motion. Then, at a random time T which verifies 
^ r{Ys)ds ~ Exp(l) ^F{T > t\Ys, s < t) = exp r{Y,)ds^ 

it jumps following F. Indeed Xt = F{Xt-, 6) = F{Yt-, Q) where is a uniform variable 
on [0, 1]. Then, the process evolves like the diffusion and we repeat these steps again and again. 
This process is an hybrid process 151 and when a = it is a piecewise deterministic Markov 
process (PDMP) lfT4ll . 

Some subclass of this process was studied. For instance, the articles |[3l [T0l[T5l[T7l[T8ll26ll25l[32l 
study a version which represent the modeling of the famous Transmission Control Protocol 
(TCP) used for data transmission over the Internet. And in ||T]|2l[T3l|24l[3Tl, some others exam- 
ples represent the growth of some biological content of the cell (nutriments, parasites,...) which 
is shared randomly in the daughter cells when the cell divide. There is again a lot of application, 
like reliability theory and queuing llT4ll23ll25l . Some theorems about the ergodicity and others 
properties are established in these papers. But to find an explicit bound for the speed of conver- 
gence is an interesting question which it remains to establish. 

Our aim in this paper is to get quantitative estimates for the convergence to equilibrium of 
X. For some particular- coefficient, it is ergodic and admits a unique invariant law. We can find 
conditions in 1231 . when the jumps are additive, and in H] proposition 5.1], when the jumps are 
multiplicative. These two results are proved via Lyapunov techniques. Nevertheless, this pro- 
cess is, in general, irreversible since time reversed sample paths are not sample paths and it has 
infinite support. This makes Lyapunov techniques less efficient for the derivation of quantitative 
exponential ergodicity. Furthermore, another main difficulties is that entropy methods fails. In 
general, the process do not verify a Poincare or log-Sobolev inequality. 

Using different techniques (gradient estimate via Feynmann-Kac formula and non coalescent 
coupling), we arrive to prove different bounds in Wasserstein metric. These bounds are optimal 
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in term of curvature of semigroup. The proofs are inspirated by ||9j \T0\ |27l which gives similar 
result for others models or for particular cases. Recently, to find quantitative bounds, Harnack 
inequalities, functional inequalities or Wasserstein curvature, for jumps processes, has spurred 
an enormous amount of research |[3ll4ll9l[T2l[T6l[T7ll22ll23l[37l. 
The remainder is explained after the main results. 

1.2. Lipschitz contraction and Wassertsein curvature. In this section, we present briefly 
some basic definition about Lipschitz contraction of Markov process. The results presented 
here hold in a very general setting. Some of these are new, even if there is some close theorem, 
and some others are not new, but we hope our proof is simpler. Let {E, d) be a Polish space, 
we denote by Vd{E) the set of probability measures fj, on E such that /^g^; d{x, XQ)fi{dx) < oo 
for some (or equivalently for all) xq G E. We recall that the Wasserstein distance between two 
probability measures /^i, /i2 G T^diE) is defined by 



2, 



(3) Wrf(/Lti,/i2) = inf // d{x,y)u{dx,dy), 

!^6Marg(/ii,/i2) J J E'^ 

where Marg(^i, /i2) is the set of probability measures on E^ such that the marginal distributions 
are ni and /U2, respectively. This infimum is attained [34 j . The Kantorovich-Rubinstein duality 
041 Theorem 5.10] gives that we also have the following representation: 

(4) Wd(w,/^2)= sup gdfii- gdfi 
where Lip(d) is the set of Lipschitz function g with respect to the distance d, i.e. 

II II \g{x) - g{y)\ ^ 

ll5llLip(d) := sup < oo, 

and hipi{d) = {g £ hip{d) \ ||5'||Lip(d) < 1}- If the semigroup is smooth enough, we have 
5xPt G 'Pd(M+), for all X G and t > 0. Hence, the semigroup is well defined on Lip{d). 
The Wasserstein curvature of {Xt)t>o with respect to a given distance d is the optimal (largest) 
constant p in the following contraction inequality: 

(5) l|^t||Lip(d)^Lip(d) < e-f\ t > 0. 

Here ||-Pt ||Lip(d)-^Lip(d) denotes the supremum of ||-Pt/||Lip(d) when / runs over Lipi((i). It is 
actually equivalent to the property that 

This notion of curvature was introduced by Joulin Il22ll and Ollivier |[29ll and is connected to 
the notion of Ricci curvature of Riemannian manifold. This definition of curvature is relatively 
close to the notion of Wasserstein spectral gap (see for instance, |[T9l ) : A semigroup possesses 
a Wasserstein spectral gap if there exist A > and C < such that 

m{S.,Pt,6yPt) <Ce-^'d{x,y). 

The presence or not of the constant C may be important for concentration inequalities. If the 
optimal constant is positive, and the process has a stationary probability measure vr then the 
semigroup converges exponentially fast in Wasserstein distance Wd to the stationary distribution 
(it is enough to take /i2 = ). We can notice that an exponential decay is possible even though 
the Wasserstein curvature is null: 
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Theorem 1.1 (Exponential decay when p = 0). Let X be a (general) Markov process such that 
its Wasserstein curvature defined at Q is null (or non-negative). If there exists to > such that 

Kto - sup — < i. 



Then there exists k > 0, such that for all xq, yo £ 



Furthermore, k can be chosen as 

HKto) 



K 



to 

The Wasserstein curvature is a local characteristic. It corresponds to the worst exponential 
decay. In this paper we use two different distances. When d{x, y) = \x — y\ for all x,y £ E, 
Wd is the standard Wasserstein distance (also called Kantorovich distance). So, in this case, we 
denote it by W and d = | • |. When C M and Fi, F2 are the cumulative functions of /ii and 
fi2, we have 

W(/ii,/i2)= / \Fi{x)- F2{x)\dx=\\Fi-F2\\Ll. 

Je 

The convergence in Wasserstein distance is equivalently to the convergence in law and the con- 
vergence of the first moment. 

When d{x, y) = '^{x=y} for all x,y e E, Wd is the standard total variation distance. So, in this 
case, we denote it by dxv- The equation ^ becomes 



djwilJ,i,P2) = ^ sup / fdi^ii - [ fdfi2- 

^ II /■Il_<1 JE JE 



ll/l|oo<l' 

If [i\ , /i2 have density n\ , 712 with respect to the Lebesgue measure, we have 

1 /■ 1 

diy{[i\,[i2) = 1: \ \n\{x) - n2{x)\dx = -\nx- n2\L-^. 
I JE 2, 

We finish this section with two results. The first one gives the existence of an invariant distribu- 
tion and the second one describes the link between Wasserstein decay and L^— spectral gap. 

Theorem 1.2 (Existence and uniqueness of the stationary distribution ). Let {Pt)t>o be a Markov 
semigroup, if its Wasserstein curvature p is positive then there is a unique invariant measure vr. 
Furthermore, 

Vi/ G Vd{E),Wd {pPt, vr) < e-P'Wd {u, vr) 
We can compare this theorem (and its proof) with ifTTl Theorem 5.23]. 

Theorem 1.3 (Wassertsein contraction implies L^— spectral gap for reversible semigroup). Let 

iPt)t>o be a semigroup of a Markov process. Assume: 

• Its Wasserstein curvature, p, is positive. 

• The semigroup is reversible. 

• Lipi{d) n L°^(7r) n L'^{tt) is dense in L'^{tt). 

Then there is a L^— spectral gap. Indeed, 

Var^iPtf) < e~^P'VarM)- 
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Ptf - j fdn 




f- 1 fdTT 













Or equivalently, 



We can see, in the proof, the importance of the reversibility. Nevertheless, for the most 
part of our examples, this condition is never verified. This theorem is just the continuous time 
adaptation of |[T9l Proposition 2.8]. We can compare it (and its proof) with [11, Theorem 5.23] 
and Il36l Theorem 2.1 (2)] for continuous time version and ||29l Proposition 30]. Notice that 
the two last theorems are also true when there is a Wasserstein spectral gap instead a positive 
curvature. Let us end this section by advertising on a parallel and independent work |f33l which 
give some related results. 

1.3. Main results : Quantitative bounds of jump-diffusions. In this section, we first recall 
some properties of the process that we study, then we state our main results. Let Ptf{x) = 
E[f{Xt)\Xo = x] be the semigroup generated by O. By the Ito-Dynkin, it means that 



(6) 



Ptf{x) = fix) + f P,Cf{x)ds = fix) + f CPsf{x)ds 
Jo Jo 



for all / smooth enough (see or llT4l for the domain of the generator). We can describe the 
evolution of X in terms of stochastic differential equations (SDE). Let (i?s)s>o be a standard 
Brownian motion, and let Q{ds, du, dO) be a Poisson point measure on M+ x M_|_ x [0, 1], of 
intensity ds du dO, independent from the Brownian motion. Here, ds, du, dO are the Lebesgue 
measures on M4., R+ and [0, 1]. Then we have. 



Xt = Xo+ f g{Xs)ds + /* V2a{Xs)dBs 
Jo Jo 

l{«<r(x,_)} {F{X,.,e) - Xs^) Q{ds, du, dO). 



(V) 



IE x[0,l] 

We are interested by the long time behavior and not by the existence. Hence, in all the paper, 
the coefficient are supposed to be smooth. Now, we are able to express our main results. 

Theorem 1.4 (Wasserstein curvature for jump-diffusion which conserves the monotonicity). 
Assume ^ ^ 

j d.^F{x,e)de ^Qandr'{x) [ x - F{x,e)de < 0. 
Jo JO 

-g'{x) + r{x 



If 



inf 

x>0 



Fix, 



p>0. 



d^F{x,e)dej - r'{x) 

Then the contraction inequality ^ is satisfied with the optimal constant p. Indeed, we have 

W{pPuuPt)<e-P'W{p,v) 
for all measure p and v which have a first moment. 

In section ??, we will see that, under this condition, the process is stochastically monotone. 
Indeed, the semigroup preserves the monotonicity. It is easy to see that the conserve is true. In 
particular this theorem can be used when 

• r is decreasing. 

• Fix, 9) < X, for almost all 6 G [0, 1]. 
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• X ^ F{x, 9) is increasing. 

For instance, the jumps can be multiplicative (indeed F{x, 9) = ip{9)x, where ip < 1), additive 
(indeed F{x, 9) = x — ip{9), where 99 > 0), or a combination of these two types. We can, of 
course, also suppose that r is increasing, and, for almost all 9 € [0, 1], F{x, 9) > x. 

In the expression of p, we can see the interplay between the drift parameter (indeed g) and 
the jump mechanism (indeed r and F). We also see that the curvature do not depend to the 
diffusive term a. We also give an other criterion to find the Wasserstein curvature, which is just 
a generalisation of the result of ifTOl theorem 2.3]: 

Theorem 1.5 (Wasserstein curvature when r is affine). If the following assumptions hold 

• E cR+. 

• F{x, 9) < X, for almost all 9 G [0, 1]. 

• X ^ F{x, 9) is increasing. 

• d,Fix,9) < 1. 

• r(x) = ax + b, for all x >0, where a,b > 
Then the contraction inequality (O verifies 

P> -supg'(z) + 6(l-sup / d^F{z,9)d9] . 

This bound looks like the bound of theorem [L4l (because E C M+ and F{x, 9) < x impUes 
inixe E x—F{x, 9) = 0). It will be interesting to prove that this bound is optimal and right for all 
r. These two theorems give interesting estimates in Wasserstein distance. It is also interesting to 
prove some bounds in total variation distance, but it is more difficult. For instance, for the TCP 
windows size process, between the jumps the trajectories are parallel, thus it seems impossible 
to intersect it (see section |4|. A first approach is given in |l3|. 

It will be interesting to find the same type of result when we add a discrete component to our 
process. It will be generated by 

Lf{x, n) = a{x, n)f"{x) + g{x, n)f{x) + r(x, n) (^j^ f {F{x, n, 9)) d9 - f{x, n)^ . 

In this case, X = {Z, N), and Z may be continuous (see |[T6l . for instance). It is also interesting 
to investigate the multidimensional case. The proof of theorem [L4l is not adaptable in these two 
cases. 

The papers is organised as follow: The next section is devoted to the proof of theorem about 
Lipschitz contraction. The proof of theorem [L4l and theorem [T3] are in section |3] Many exam- 
ples and applications are given in the last section. 

2. General properties about Lipschitz contraction 

In this section we prove the theorems of section 11.21 

2.1. Exponential decay when p = (proof of theorem ll.ll l. Let X be a Markov process such 
that its Wasserstein curvature defined at ^ is null (or non-negative). Assume there exists to > 
such that 

Kto - sup — < i. 
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We will prove that there exists k > 0, such that for all xq , yo G E, 
Let uj be defined, for all x y,by 



uj{t,x,y) 



and 



Wd{£{Xt\Xo = x),£iXt\Xo = y)) 
\x-y\ 

u){t) = sup uj(t, X, y). 



Let (y, Z) be any coupling of the law of X. Indeed Y = Z = X but Y and Z are not necessary 
independent. The Markov property gives 



y, z) = sup E [/(yt+,) \Y^ = y]-^ [f{Zt+s) \ Zo = z] 

/GLipi 



= sup E E f{Yt) \Yo = Y, - E 
/eLipi ^ ^ ^ 

<E[LD(t)d(n,Z,) |yo = y,^o = ^] 



/(Zt) I Zo = \Yo = y,Zo = z 



where {Y, Z) is a coupling of the law of X. Taking the infimum on all coupling we deduce, 

u{t + s,x,y) < uj{t)uj{s,x,y) u}{t + s) < u}{t)u]{s). 
Now, by the hypothesis, the curvature is non-negative, thus 

yt > 0, uj{t) < 1. 

Hence, ui is decreasing. But, by the assumptions, there exists to > such that co{tQ) < 1, and so 

Vt > to, ^(i) < 1- 
Finally, for all t > Iq, there exists n G N such that t > nto, and then 
Wd{6xPt,6yPt) < co{t) m{5.,6y) < uj{n^) Wd{5x,6y) 



< u{tor mSx, 5y) < exp [t l^^MMZ ^ Wd{5x. 5y). 



2.2. Existence criterion of the invariant distribution (proof of theorem 11.21 1. Let X be a 

Markov process and {Pt)t>o its semigroup. First, let us assume that 

. Wdil^Pt, uPt) < e-P'Wd{fi, y) for all £ Vd{E). 

• C = sup^>o VVrf((5a;, ^xPt) < +00 foT somc X e E. 
We will prove that gives the existence of an invariant measure. The proof which follows is 
simple and comes from a discussion with D.Villemonais. He has developed the same type of 
criterion, for the existence of Quasi-stationary distribution, with the help of P. Del moral in a 
work still to be completed (see i35i chapter 5]). Consider the sequence {6xPt)t>o, it verifies 

yt,s> 0, Wd{6xPt+s,6xPt) < e-P''Wd{SxPs,6^) 

< e-P'C. 
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Furthermore, the space ('Prf(i?), W^) is known to be complete lITTl theorem 5.4] thus the se- 
quence {dxPt)t>o converges to a measure vr G Vd{E). Then 

lim Wd{vPt+s,r]Pt) = Wd{T^Ps, vr) = 0. 

Thus vr is invariant, and the Lipschitz contraction gives the uniqueness and the exponential decay. 
Now, it rest to prove that sup^>o V^^di^x, ^xPt) < +oo for some x e E. We have 

sup d{Ptf{x),Ptf{y)) < e-P'd{x,y). 

ll/l|Lipi<l 

Hence Pt remains stable the space of Lipschitz function. As we have \Vd{^x,^xPt) = Ptfx{x), 
where fx{y) = d{x, y), the mapping x H> Wdi^x, ^xPt) is continuous and for all x G -E and 
e > we have 

(8) sup Wd{5x, 5xPt) <+oo. 

t6[0,e+l] 

And we have for all A; G N* and t£ [e, e + 1], 

fe-i 

m{5x,5xPkt) < Y.^di5xPjt,SxPij+i)t) 
i=i 

fe-i 

<J2e-^P*Wdi5x,SxPt) 
i=i 

< Wd{Sx,6xPt)^—^ 



The last inequalities and ^ gives the needed boundedness. 

2.3. Lipschitz contraction and L^— Spectral Gap (proof of theorem 11.31 1. The proof of this 
result is a continuous-time adaptation to [il9i Proposition 2.8 ]. Let / be a non-negative, Lips- 
chitz and bounded function such that fdir = 1. Using the reversibility, we have, 

Var^(Pt/) = l^fP2tfdTr-(^J^fd7ry 
< ll/llLip{d)W,(P2t/f^vr,7r) 

where P2tfd'K is the measure which verifies, for all smooth 99 

/ ipP2tfd-K= / / ip{x)f{y)P2t{x,dy)TT{dx) 
Je JeJe 

= v{x)f{y)P2t{y,dx)-K{dy). 
JE JE 



Because vr is reversible. Hence, 

Var,(PJ) < ||/||Lip{d)Wrf((/7r)P2t,vr) 
(9) < 
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By translation/dilatation, the last inequality holds for all Lipschitz and bounded /, and by 
density, for all / G ^^(tt). Now, let / be a measurable function such that, 

f{x)TT{dx) = and / f{xfTT{dx) = 1. 

E JE 

Applying spectral theorem, Jensen inequality and dill, we find 

/r roo 
PtfdTT = / / e~^'dEx{f)d7r 
JE Jo 

Taking the limit s — )• +oo, we conclude the proof. 

3. Exponential decay for jump-diffusions 

In this section, we prove the announced results about jumps-diffusion and add some remark 
or corollary. 

3.1. Gradient estimate (proof of theorem 11.41 ). In this section we follow the same approach 
to ||9l. Using a Feynmann-Kac semigroup, we estimate the derivative of our semigroup. So, we 
begin to prove, for any smooth enough function /, 



(10) 



Vx > 0, (PtfYix) = E 



Where y is a Markov process generated by 

Csfix) = cj\x)f"{x) + {2a{x)a'{x)+g{x)) f{x) 



+ r{x) fd,F{x,e)de{ I / d^F{x,e)f{F{x,9))de-f{x) 

Jo \Jf^ d,j:F{x,e)de Jo 



-r'(x) C x-F{x,e)de\—, i C r 

Jo X- F{x,e)de Jo Jpixfi) 

Lemma 3.1 (Gradient estimate via Feynmann-Kac formula). Assume 



f{u)du dO - f{x) 



If 



where 



d^F{x, e)de / and r'{x) [\ - F{x, e)d9 < 0. 
Jo Jo 

exp(^- J^V{Ys)ds^ \Yo = x 



yxG E,K 



< +00. 



V{x) = -g'{x) + r{x) (^1 - d^F{x, e)de^ - r'{x) j\ - F{x, 
Then ^ holds. 
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Proof. First, we have (£/)' = {Cs — V) f. Then, by the Ito-Dynkin formula Q, 

{PtfYix) = fix) + [\cpjy{x)ds 

Jo 

= f{x)+ f {Cs-V){Psfnx)ds 
Jo 



Thus t I— (Ptf)' verifies this partial differential equation: dtu = [Cs — V) u. But, it is known 
that the following semigroup verify this equation too: 



Stf{x) = E 



where Y is generated by Cs- As there is a unique solution, we get {Ptf)' = Stf ■ 



□ 



Corollary 3.2 (Propagation of monotonicity). Under the same assumption, if f is non-increasing 
then Ptf is also non-increasing. 

This corollary implies that, for all x < y, there exists a coupling such that the 

marginals are generated by C, which start from (x, y), and 

\/t > 0, Xf < X\ a.s.. 

It is why this type of formula is available only for a non general class of r. 

proof of theorem lL4\ Let / be a Lipschitz and smooth function. For any x, y > 0, we have by 
the intertwining identity (ITOl ). 

\Ptf{x) - Ptf{y)\ < sup |(Pt/y(z)||x - y\ 



z>0 

< supE 

z>0 

< supE 

z>0 



\x - y\ 



\x - y\ 



so that dividing by |x — y| and taking suprema entail the inequality: 

\\Pt\\Lip->Lip < sup E exp f- / V(Ys)ds] \Yo = z 
z>o L V JO / 

Finally, the right-hand-side of the latter inequality is nothing but ||Pt(id)||Lip, showing that the 
supremum over Lip is attained for the function id : x ^ x. It achieves the proof. □ 

Remark 3.3 (/i— transform and first eigenvalue). Assume that Cs — V have a first eigenvalue 
such that its eigenvector ^ is positive. Using an h— transform with h = e~^*^, we get for any 
Lipschitz function f, 



E 



e"-^V(2;)E 



\nzt)i 


< e"^*E 




[HZt)\ 




mzt)\ 



Then, ifip is smooth enough, the Wasserstein spectral gap is A. See section 11 for an application. 
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Theorem 3.4 (Exponential decay when the curvature is null). Under the assumption of the 
previous lemma and ifY is irreducible, V > 0, and there exist O = (a, h) and e > such that 
is compact and 

Vx € 0,V{x) > e. 
Then there exist to ^ and k > such that 

Vt > to, uPt) < e-"* v) 

for any n,^ £ ^11(1^+). 

a and b may be infinite if there is not in the same time. Notice that, if there exist x £ E, such 
that V{x) = then the curvature is null. 



Proof. The proof is adapted to 127] section 5.1]. Let 



D{t,x) = 



exp 



V(Ys)ds 



and D{t) — sup^>Q D{t, x). It is easy to see that for all x G E, D{-, x) and D are decreasing. 
Furthermore, D{t, x) = 1 if and only if, starting from x, Ys G almost surely, for almost all 
s < t. Then, as Y is irreducible, we have 

(11) Vx > 0, Vt > 0,L>(t,x) < 1. 

If X G O, then it is also right for t = 0. We begin to prove that there exist to such that 

Vt > to, D{t) < 1. 

Then, we will prove the assertion of the theorem. For x G O, set r = inf{t > | G O^}. We 
have. 



D{t,x)=E^ 



lr<t exp 
lr<t exp 



V{Ys)ds 
V{Y,)ds 



lr>t exp 



ViYs)ds 



And 



lr<t exp 



V{Ys)ds 



< E. 



< E. 



< E^ 



+ e 



-et 



exp ^ V{Y,)ds^ I J'r 



lr<te max (D(t-r, c)) 

ce{a,b} 

lr<V2e"^* max {D{t - r, c)) 

ce{a,b} 



h./2<T<te " max (D{t-T,c)) 

ce{a,b} 



< max (Z)(t-r,c)) +e-^*/2 

c£{a,b} 
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where T is the fihration associates to Y . Thus, 



supL'(i,x)< max (D(t-r, c)) + e ^ +e 



-et/2 



Then we deduce that limsupj_^^oo sup^^go x) < 1 and then the existence of to such that 
for all t > Iq, 

sup D{t, x) < 1. 
xeo 

When X G O'^, we see that he Feynmann-Kac semigroup is continuous because Y verifies the 
Feller property, then, the supremum of sup^,g0c D{t,x) is attained. And finally, by (fTTT) . we 
have D{t) < 1 for all t > to- Now, we can use theorem [T7T] to conclude. We can also use the 
argument to ETl . That is, the Markov property ensures that 



D{t + s,x) =E^ 
<E, 



exp 



exp 



V{Y^)duj Ey, 
V{Yu)du) D{s) 



exp 



V{Yu)du 
D{t,x)D{s) 



< D{t)D{s) 



and then D{t + s) < D(t)D{s). For all t > to. there exists n G N such that t > nto- Hence, we 
have 



WifiPt,i^Pt) < D{t) W{ii,v) < D{nto) 



< D(to)" < exp t 



to 



□ 

3.2. Contraction by coupling (Proof of tlieorem 11.51 ). As in fTp], we make a coupling which 
favour the simultaneous jumps. Because, by the properties of F, during a simultaneous jump, 
the distance is decreasing. More precisely, this couple is generated by 

Gf{x, y) = a{x)dxxf{x, y) + a{y)dyyf{x, y) + g{x)d^f{x, y) + g{y)dyf{x, y) 
+ r{y) f' fiF{x,e),F{y,e))de 





+ (r(x) - r{y)) / / {F{x, e),y) d9 - r{x)f{x, y) 
Jo 

if X > y, and 

Gf{x, y) = a{x)dxxf{x, y) + a{y)dyyf{x, y) + g{x)d^f{x, y) + g{y)dyf{x, y) 
+ r{x) r f{F{x,e),F{y,e))de 





+ {r{y) - r(x)) / / (x, F{y, 0)) dO - r{y)f{x, y) 
Jo 

in the other case. The first line describe the space motion between the jumps and the second line 
the jump mechanism. The dynamics of the couple of components is as follows: 
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(1) After an "appropriate" time which depends of the two trajectories (or only one in some 
cases). The upper coordinate jumps. 

(2) Simultaneously, the other one "tosses an appropriate coin" whose probability of success 
depends on the positions on the two components to decide whether or not it jumps too. 

(3) In the case of joint jumps, both components use the same uniform variable. 

(4) Then, we repeat these three first steps again and again... 

Let a be defined by a(j. = E[|X( — Yt\], where {X, Y) is generated by G and starting 
from {x,y). It verifies j^)(0) = Gf{x,y) where / : {u,v) ^ \u — v\. Assume x > y, we 
have. 



a 



[,,y){0)=9{x)-g{y) + {ay + b) / {F{x,6) - F{y,9))de 







1 



+ a{x — y) I \F{x, 9) — y\d9 — {ax + b){x — y) 
Jo 

<{x- y)[supg'{z) + {ay + b) sup / d^F{z,9)d9 

+ a \F{x, 9) -y\d9- {ax + b)\ 
Jo 

< {x-y)[snpg'{z) + blsup d^F{z,9)d9 - l] 
z>0 \zeEJo ) 

+ a f' If{x, e)-x + y{l- sup f' d,F{z, 9)d9)] lF(x,e)>y 
Jo \ z&eJo J 

+ (^y-x)+y{l- sup d^F{z, 9)d9)^ lF^^^e)<yd9] 
<{x-y) sup g'{z) + b (sup / dccF{z, 9)d9 - 1 

\_z>0 \zeEJO 

Then, by the Markov property, 

1 ) ) 



a\^y){t) < (sup 5' (z) + 6 (sup f^d^F{z,9)d9 

^ \z>0 \zGEJ0 

^Oi{,x,y){t) < exp ^supfii' + b ^sup^ d.xF{z, 9)d( 



l\t\\x-y\. 



The end of the proof is standard. 



4. Examples and Applications 

We begin by the easier example. This example show the interplay between the growth and 
the jump rate. After, we consider five examples. The first one represents a process, which grows 
Unearly, and jumps additively. Then, the three next examples are different growth-fragmentation 
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processes. We called growth-fragmentation process, a process generated by 

(12) Cf{x) = a\x)f"{x) + g{x)f'{x) + r{x) (^J^' f{hx)n{dh) - /(x)) . 

The state space can be -E = M or £' = M+. And we finish by an example of structured- 
population. 

4.1. Straightforward consequences. 

4. 1. 1. Jumps independent of the current state. Assume 9) = F{6) for almost all 6 G [0, 1]. 
It means that, when a jump occurs, the new position of the process do not depend of the previous 
position. This process forgets its origin just after its first jump. We can find an expression of 
the invariant distribution in [^51 . In this case, the theorem [L4l is valid when r increases before 
/q^ F{6)d6 and decreases after. And we have that the curvature p verifies 

p > inf r(z) — supg'{z). 

4.1.2. Feller diffusion with multiplicative jumps : The Bansaye-Tran process. Let us consider 
the process studied in Q. It evolves like a Feller diffusion between the jumps. Indeed, it moves 
following the E.D.S.: 



Vi > 0, dYt = gYtdt + y^2aYtdBt. 

Where g are o" are two positive numbers. And, when it jumps from x, this new state is Hx 
parasite, where H ^ T-L, and is a measure supported in [0, 1]. This process models the rate 
of parasite in a cell population. The parasites grow in each cell and, sometimes, the cells divide. 
These two phenomena do not unwind in the same time scale. The parasite are bom and die faster 
than the cell divide. Hence, the rate of parasite is modelled by a Feller diffusion. This one can 
be understood as the limit of birth and death process. The jumps models the division of cells. In 
this setting, we have 

Corollary 4.1 (Exponentially decreasing to 0). Ifr is decreasing or ajfine and 

g< inf r(x) / h'H{dh) 

then 

E[|Xt|] < exp (^-t (^^[r{x) hn{dh) - g^^ E[Xo]. 

Proof. Under theorem [L4l and theorem [T31 the Wasserstein curvature is positive. Furthermore, 
5o is invai-iant and Wid^Pt, Sq) = Ptld{x) = E[Xt \ Xq = x]. □ 

The results of Q are better even if this corollary gives a bound for the L^-convergence. To 
compare, for instance, EJ proposition 3.1]) tells us 

Theorem 4.2 (Extinction criterion when r is constant). We have the following duality, 
(i) Ifg < -r Jo log{h)n{dh), then F{3t > 0, = 0) = 1. 
Moreover if g < —r log{h)'H{dh), 

3a > 0, Vxo > 0, 3c> 0, P(Xj > | Xq = xq) < ce""* {t>0). 
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(ii) Ifg > -rj^ \og{h)n{dh), then P(Vt > 0, > 0) > 0. 
Furthermore, for every 0<a<g + r/ log{h)Ti (dh), 



lim e-'^'Yt = oo = {Vt, Xt > 0} a.s. 

f— >+co J 

When r is constant and g < rn, we have g < — r Jq \og{h)T-L{dh) . Thus X converges almost 
surely to 0, and 

nXtW < e-('^'^-^)E[|Xo|]. 
More precisely, a rapidly calculation (which is given in the more general case hereafter) gives 

E[Xt] =e(e"^^)*E[Xo]. 

So, if 

/ - \og{h)n{dh) > i > K = 1 - / hn{dh) 

Jo r Jo 



then 



lim Xt = a.s. but lim E[Xt] = +oo. 



These two convergences explain why our criterion is not applicable and why the convergence is 
slow in this case. 

4.1.3. Speed of convergence for branching measure-valued processes. Consider a model of 
structured population. We observe a content 1" of a cell which evolves following the stochastic 
differential equation Q and after an exponential time T, the cell dies and is replaced by a ran- 
dom number of offspring K. The content of the mother is shared in each new cells. Indeed, the 
offspring j < K have the content F{Yt, &) . Consider the measure defined, for all t > 0, by 

u&Vt 

where X^ is the content of the cell u at time t, and Vt the set of individual alive. It was proved 
in Q and HSl that 



E 



f{x)Zt{ds] 



nf{Xt)]. 



Where X is generated by ([T]), with biased parameter. Hence, we deduce the long time behavior, 
and the contraction properties of the mean measure. A same formula holds when r is not constant 
(see |[T3l ). It will be interesting to capture the speed of convergence to Z instead of E[Z]. A first 
approach is given in |[T3l theorem 1.2]. 



4.1.4. Stress release, repairable system and workload models. Assume that £^ = M, and let us 
define 

Cif{x) = gf'ix) + r(x) ['fix- Gix, e)) - fix)d9, 

Jo 

where g is a positive number and G is a non-negative and measurable function. The process 
generated by £i was studied in |[23l . This class represents a generic model of applied probability 
and is of importance in earthquake modeUng, reliability theory and queuing. We have 
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Corollary 4.3 (Wasserstein curvature when r is decreasing). If r is decreasing, x G{x, •) is, 
almost surely, increasing, we have 

p> inf (r{z) d^G{x,e)de] -supr'(z) G{x,e)de. 
^eR \ Jo ) zm JO 

In particular, if G do not depend to x, then we have 

p = -supr'{z) G{e)de. 
^gk Jo 

This corollary implies a bound to the Wasserstein curvature of the virtual waiting time (also 
workload) process. Indeed, this process is generated by 

C2f{x) = -gl.>o/'(x) + r{x) / (x + G(x, 6)) - f{x)de 

Jo 

where g is a positive number and G is non-negative. If W is generated by C2, it is a lipschitz 
transformation of a process generated by £1. Indeed, we can write 

Vt >^,Wt = max(Xt,0). 

where X is generated by £1 (with an other G). 

4.2. Kolmogorov-Langevin processes. Let us consider the process which verifies the follow- 
ing S.D.E.: 

dXt = V2dBt - q'{Xs)ds 
where q is C°° and S is a standard Brownian motion. It is already known, under suitable assump- 
tion, that this process converges to the Gibbs (or Boltzmann) measure ^{du) = e"'^^'^^ du/ Z 
(where Z is a renormalizing constant). The theorem 11.41 shows that the curvature is equal to 
p = inf^gK q"{z). This result was akeady known but we can hope an exponential decay to the 
invaiiant measure when p < 0. The case where p = was studied in ETl . Our main example 
of interest is 

q{x) = x'^ — x^. 

We begin by a preliminary lemma 

Lemma 4.4 (Existence of eigenelements of a certain operator). Let us define H the closure on 
of the operator defined by 

Hf = -f"{x) + f^ 
for all smooth enough f.Ifq is C°° and 

lim q'{x) = +00, 

|x|^+oo 

there exist a unique positive function 99 G L^nC°° and a real number A > such that Hip = Xip. 
Proof. It is know from |l6j theorem 3.1 p. 57] that if 

iim = +00 

\x\^+oo 2 

then H has a discrete spectrum. Furthermore, there exist a unique positive eigenvector. It 
corresponds to the smaller eigenvalue A > and we denote it by ip (see ||6j chapter 2]). The 
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regularity of ip comes from the regularity of q and the uniqueness of the solution of the linear 
equation. □ 

Theorem 4.5 (Rates of convergence). Under the same assumption of the previous lemma, we 
have for all x,y gM., 



W {S^Pt, 6yPtf) < e-^* sup e'?(^)/2 f^^ ^ Ca/?*^ 



ip{z) 



\x - y\ 



where Ci, C2 are two positives constants (Ci is explicit) and (3 G (0, 1). 

Proof. Let {Pt)t>o be the semigroup of X. As can be view in the proof of theorem [L4l we have, 
for all / smooth enough, 

= Stf 

where {St)t>o is a Feynamnn Kac semigroup generated by L defined by, 

Lf = fix) - q'{x)nx) - q"fix). 
Let {Qt)t>o be defined for all smooth enough / by 



We have 



Where 



dtQtU) = X—e-'^/^Stif^e'^/^) + —e~'^/^St{L{f^e'^/^)) 

ip If 

(f if 
= ^—e-'i/^St{e'i/\Gf) 
= QtiGf). 



Gf = f" + 2^/'. 



This relation can be understood as 



where 1" is a Kolmogorov-Langevin process generated by G and starting from x. Hence, 



\dM{x)\ < e 



■ inf 



But 



G- 



-if 



So, by ll28l theorem 6.1] with V = we have 



X-q" 

f 
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for all x G M, where /3 G (0, 1), > and 7r{dx) = ef ^'^'^'^dx/Z = ip{x)^dx/Z {Z is a 
renormalizing constant). Finally, 



□ 



Example 4.6 (Omstein-Uhlenbeck). If q{x) = f/ien f/je theorem is not available. But we 

see that ip = lis an eigenvector ofH with respect to the eigenvalue /x/2. It is the smaller eigen- 
value such that the eigenvector is positive. With a same computation to the previous theorem, 
we can prove 

W{6,Pt,6yPt)<e-^''/^\x-y\. 
We can compare this result with the one of theorem U .4\ 

W{6^Pt,6yPt)<e-^''\x-y\. 

In fact, it easy to see that it is an equality. 

This example prove the non-optimality to the previous theorem. But as can be viewed by the 
following example, this techniques can be useful for more complicated example. 

Example 4.7 (An example where q is not convex). Let us consider 

q{x) = — x'^. 

Hence q"{x)/2 = — 1 and Hip{x) = Xip{x), when 

X = — 1 and (p{x) = exp (^—x'^VQ/2^ . 
So, as vr, the invariant measure, integrates 1 /ip we have, 

>V(/iPt,7r) <C^,,e-(^-i)* 
where ji is a measure which integrate Xjip and C^,7r ci constant which depend to ^ and vr. 
4.3. TCP Windows size process. 

4.3.1. An ideal model for the TCP windows size. Now, we study some model which represent 
the TCP congestion. The first one was developed in |[30ll . Before a jump, the process evolves as 

dtXt 



The jumps arise following a Poisson process of intensity r with points (T/^)fc>o- And the 
jumps are multiplicative, indeed, Xt^ = IiXt^^. The parameters verify 

< h < 1, m > -1, g > and r > 0. 

In QUI , the proof of the ergodicity and an explicit formula for the invariant measure is proved: 



Theorem 4.8 (Convergence to equilibrium). For any initial condition, X converges in law to a 
measure vr which has a density defined by 



X ^ ^ afc(h) exp —x^^^ 



k>0 



;(m + 1) 
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where 

"'^^ ^ ~ ut=ii^ - ni=i(i - hi("^+i)) ■ 

By theorem [L4l if m > 0, the Wasserstein curvature is p = r(l — h). 

Corollary 4.9 (Wasserstein exponential ergodicity). Assume m > , for any probability mea- 
sure G 7^|.| (M+) (indeed, which has a first moment), we have 

4.3.2. A more realistic model. Now, we consider that the jumps depend of the position and the 
growth is linear. Let {Xt)t>o be the Markov process generated by 



(13) 



Vx > 0, Cf{x) = fix) + rix) f{hx)n{dh) - f{x)^ 



To fix an idea, r{x) = x is an interesting value |l3][l5l[T8l, but, the case where r is constant, is 
often used for its simplicity |[26l l30l . The invariant measure is known when r(x) = rx"(it is 
explain in the following). We have: 

Corollary 4.10 (Wasserstein curvature). Ifr is decreasing, we have that p, defined at (|5]), verifies 

p = [I — / liHidh) ) inf (r(x) — xr' (x)) . 
\ Jo J ^>o 

If otherwise, r(x) = ax + b, we have 

p> (^^- hn{dh)^ h=[l- hn{dh)^ inf {r{x) - xr'{x)) . 

In particular, when r{x) = x, we have p > 0. When the process is near than 0, the jump rate 
is null. Hence the process evolves like t i— > Xq + t which have a null curvature. So we see that 
p = 0. A recent work ||3l prove that, nevertheless the curvature is null, this process converges 
exponentially rapidly to its invariant distribution. 

4.3.3. The embedded chain. We always consider a process generated by ([T3] ). Let iXn)n>o be 
the embedded chain. It is defined by, 

(14) Xn = Xt„ where r„ = inf{i > T„_i | Xt+ / Xt-} for n > 1 and Tq = 0. 

This Markov chain is easier to study than the continuous time process. For instance, for r{x) = 
ax", it is easy to see that 

R{Xt„^,) = R{Hn{R~\En + R{XtJ))) = K^\En + R{XtJ) 

where R is the antiderivative of r. Then, it is ergodic and its limit Xoo verifies 

r{XtJ = h:+\e^^ + r{x^)). 

Now, using ifTSl proposition 5], we deduce that Xoo have a density given by 

^ ^ t TT " ^ -h-("+i)("+i)a(a+l)-ix("+i) 

iln>l\^ " J n>Ok=l 
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Now, applying [14, theorem 34.31], we can deduce the invariant law of the continuous time 
process. This result generalises ifTOl Theorem 2.1] but it is already known via others techniques 
|[32l . For this Markov chain, we arrive to bounded all Wasserstein distance. Indeed, recall that 
for every p > 1, the W^^^ Wasserstein distance between two laws /ii and fi2 on E with finite p"^ 
moment is defined by 

W(^)(w,/^2) = ^^ inf (E[|X-yn)^/^ 
where the infimum runs over all coupling of fii and ij,2. 

Theorem 4.11 (Wasserstein exponential ergodicity for the embedded chain). Assume that C{Xq) 
and C{Yo) have finite p'^ moment for some real p >1 and r is increasing. Let X and Y be the 
embedded chains of X and Y. Then, for any n > 0, with a random variable H ^ %, 

>v(p)(£(x„),£(y„)) < E(Fn"/^'Wp(£(Xo),/:(yo)). 

In particular, if it is the invariant law of X then 

W(^')(£(X„),7r) < E(Ff)-/f>Vp(£(Xo),^). 

Proof. It is sufficient to provide a good coupling. Let x > and y > be two non-negative 
real numbers, and let (-En)„>i and {Hn)n>i be two independent sequences of i.i.d. random 
variables with respective laws the exponential law of unit mean and the law T-L. Let X and Y be 
the discrete time Markov chains on [0, oo) defined by 

Xq = X and Xn+i = Hn+iR~^{R{Xn) + En+i) for any n > 
Yo = y and Y^+i = Hn+iR-HR{Yn) + E^+i) for any n > 0. 

the law of X (respectively Y) is the law of the embedded chain of a process generated by 
L and starting from x (respectively y). Now, let a be a non-negative number, if tpa : x ^ 

R-~^{a + R{x)) then 

(15) ^'a{x) = =^ l^aix) - iPa{y)\ < \x - y\. 

r[a + X) 

And we get 

Vp > 1, E[|l„+i - = E[F^+,|<^s„_,,(X„) - ipE,.^,iYr,)n 

< m'n+l\^n - Yn^ = E[HP_^,n\Xn " ^k^] 

A straightforward recuiTcnce leads to 

E[\Xn-Yn\P]<E[H^r\x-y\P. 

This gives the desired inequality when the initial laws are Dirac masses. The general case follows 
by integrating this inequality with respect to couplings of the initial laws. □ 

The fact that, between the jumps, the growth is linear is fundamental in this proof. It is why, 
we do not generalise it. So, we arrive to prove a convergence criterion for the continuous time 
process when r decreases and for the embedded chain when r increases. We are not able to 
associate these two approaches. 
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4.4. Integral of Levy processes with respect to a Brownian motion. Before to talk about 
integral of Levy process, we give an example of growth-fragmentation processes, where some 
calculates are explicit. 

4.4. L Generator which conserves the polynomial functions. Consider the process generated by 
([12] ) with the following parameters: 

g{x) = 50 + 9ix, a{x) = ctq + aix + cj2X^. 

And r is constant. In this case, the generator conserves the polynomial functions and their 
degree. We then deduce. 

Theorem 4.12 (Moments estimates). Let t > and x > 0, we have 
E[Xt I Xo = x] = (x + -^^) e(^--)* - 



gi-VKlJ 91- TK 

and 

Vn > 2, E[Xf I Xo = x] = + E 4"^e('^2fc(fc-i)+9ifc+r--.)t 

k=l 

where Kn = I — J h'^H{dh) and for all k G {1, .., n — 1} we have 

= Mn - 1)4"-') + (ain(n - 1) + 50^)4""'^] 

a2n[n - 1) + gin - rKn 

4"^ = 7 T(^on(n - 1)4"-') + (ain(n - 1) + <7on)4"-')) 

1 - {a2n{n - 1) + gin - rKn) 



n-l 

4") = x-E4"^- 

k=0 

Proof Let t > 0, n E N and a„(t) = E[Xf ], we get, 

ao{t) = 1, a'i{t) = go + {gi - rKi)ai(i) 

and, 

"n(*) = <Ton(n- l)a„_2W + i(Tin{n - I) + gon)an-i{t) + {a2n{n- I) + gin- rKn)an{t). 
Then the result follows by a straightforward reccurence. □ 

Remark 4.13 (Convergence to the equilibrium). This result gives the convergence of the mo- 
ments when gi — rn < 0. In this case, the n'^ moment converges to Aq^^ which is easily 
calculable. When this sequence verifies the Carleman criterion, that is 

^(4"VV2n^+^ 
n>0 

We have the convergence, in law, of X to a measure with moments defined by the sequence Aq. 

Example 4.14 (A first link with the Levy process). The If a^i = a2 = and go = 0, then X 
is the exponential of a Levy process L. The condition gi — rn < implies that X converges 
to almost surely. The Wasserstein decay of the theorem \L4\ gives a speed of convergence in 
L^-norm. 



22 



BERTRAND CLOEZ 



This example is near to the example of theorem 1421 In the following, we give an other link 
with the Levy processes and an other example which have a different behaviour. 

4.4.2. Integral of Levy processes with respect to a Brownian. Heuristically, if you take the log- 
arithm of a growth-fragmentation process, you transform the multiplicative jumps on additive 
jumps. Then, as the jump times are Poissonian, we can obtain a continuous process by renor- 
malising our process with a Levy process. Formally, let [Lt] be the Levy process defined by, 

Lt = — liu<rl ln(/i)Q((is, du, dh) = — In TT -f/^i 

Jm+ x[0,1] - \^/^^^ 

and let X the continuous process defined by Xt = Xje^*. We get. 



Lemma 4.15 (Stochastic differential equation for X and X). If X is generated by f l72l ) then, for 
any non negative number t, we get 



and, 



Xt = Xoe-^'+ f e-^^g{Xt-s)ds+ f\-^^a{Xt-s)dBs. 
Jo Jo 



Proof. All the stochastic integrals that we write are well defined as local martingales. Using 
Ito's formula with jumps EOl Theorem 5.1, p.67] (see also li2L Theorem 4.57, p.57]) and since 
X is continuous, we get: 

Xt = Xo+ f {g{X,)ds + a{Xs)dB,) 
Jo 



rt 

+ J Xs^e^'~ - Xse^''l{u<r}Qids,du,dh) 

= Xo+ f' {g{X,)ds + . 
Jo 



Then, we deduce, 

Xt = Xoe-^* + /* e^«-^* {g{Xs)ds + 
Jo 

= Xoe-^* + / e-^*-^ {g{Xs)ds + a{Xs)dBs) 
Jo 

= Xoe-^'+[ e-^'g{Xt-s)ds+ I e-^^a{Xt-s)dB,. 
Jo Jo 

□ 

This lemma is a generalisation of the relation of 12] lemma 3.2] and |[26l section 6]. In |i2J, it 
is a preliminary for the proof of theorem |42] and in |[26l . they deduce that when g is constant, 
and cr = 0, we get, 

r+oo 

Mm Xt = g e 'ds. 

The behavior of the right hand side was studied in Q and lH for general Levy processes. We 
give an other application: 
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Theorem 4.16 (Long time behavior of Integral of Levy process with respect to a Brownian 
motion). If 







where Lf = ^^=1 ^ '^^^ process independent of the Brownian motion B. Then it con- 

verges to a measure vr such that 



for all n G 2Z and 



otherwise. Furthermore, for all t > 0, 

W(£(yt),7r) <e"«*W(5o,vr), 

where k = r (^l — h'H{dh) 

Proof. In this case, we have that Y is generated by (fT2l) . where, 51 = 0, cr = 1, r is constant and 
Xq = 0. The result is an application of the results of the previous sections. □ 
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